On the quantum spectrum of isochronous potentials 
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In this paper, the quantum spectrum of isochronous potentials is investigated. Given that the 
frequency of the classical motion in such potentials is energy-independent, it is natural to expect 
their quantum spectra to be equispaced. However, as it has already been shown in some specific 
examples, this property is not always true. To gain some general insight into this problem, a WKB 
analysis of the spectrum, valid for any analytic potential, is performed and the first semiclassical 
corrections to its regular spacing are calculated. We illustrate the results on the two-parameter fam- 
ily of isochronous potentials derived in [lj, which includes the harmonic oscillator, the asymmetric 
parabolic well, the radial harmonic oscillator and Urabe's potential as special limiting cases. In addi- 
tion, some new analytical expressions for families of isochronous potentials and their corresponding 
spectra are derived by means of the above-mentioned method. 

PACS numbers: 45.05.+x; 03.65.Sq 
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I. INTRODUCTION 

Several studies have already been devoted to the in- 
teresting problem of the relation between classical and 
quantum "generalised harmonicity" . In classical me- 
chanics, the motion of an oscillator in the parabolic 
well V(x) = uj 2 x 2 /2 is known to possess an energy- 
independent frequency ui. However, this harmonicity is 
not specific to the latter and can be achieved by designing 
adequate potentials called isochronous. Their construc- 
tion simply amounts to shearing the parabolic well in 
such a way that, at fixed energy, the distance between 
two turning points is preserved (see for example 0,0). 

In quantum mechanics, the energy levels of a parabolic 
well are regularly spaced by a quantity hui. Thus, a 
generalised harmonicity is naturally defined by an eq- 
uispaced spectrum. Similar to the classical case, it is 
possible to construct potentials, essentially different from 
the parabolic well, whose spectrum is exactly harmonic. 
This can be done by applying a supersymmetric trans- 
form followed by a Darboux transform to the harmonic 
well which lead to families of potentials called isospectral 
partners of the harmonic oscillator (for a review, see for 
example Q). A famous example of such a potential has 
been derived by Abraham and Moses [j| . 

A link between these classical and quantum transfor- 
mations has been established by Eleonskii et al. @ . They 
show that the classical limit of the isospectral transfor- 
mation is precisely the isochronism preserving shear pre- 
viously mentioned. The next question is whether the two 
classes of classical and quantum generalised harmonic po- 
tentials are the same. Surprisingly enough, they are not. 
Both aspects of this problem have been investigated by 
several authors who demonstrated through specific ex- 
amples that neither isochronism generates an equispaced 
spectrum nor does a regularly spaced spectrum stem from 
an isochronous potential in general 0, B S Dll • 

Despite this inequivalence, it is worth noting that the 
semiclassical Einstein-Brillouin-Keller (EBK) quantisa- 



tion of an isochronous potential leads to an equispaced 
EBK-spectrum. Indeed, the semiclassical quantisation 
rule is Iq(E) = K(n + i), n E N where the action Iq{E) 
is related to the frequency lo(E) and the energy E by 
dlo/dE — u(E)^ 1 . As the frequency of an isochronous 
potential is energy-independent, E = loIq and EBK lev- 
els are given by E„ BK = (n + h)Huj. Thus, their spacing 
is constant, £™k „ B ebk = hjJ _ 

Consequently, as the usual semiclassical quantisa- 
tion fails to distinguish between the spectrum of an 
isochronous potential and a purely harmonic spectrum, 
studies to date have resorted either to the study of 
isochronous potentials whose quantum spectrum (or 
quantisation condition) is exactly known or to numer- 
ics. To date, the asymmetric parabolic (also called split- 
harmonic) well |38| is the only isochronous potential 
whose exact quantisation condition proves its quantum 
levels to be non strictly equidistant Q, . Another class 
of potentials denned on a finite interval of the real axis, 
and first derived by Urabe Hp . have been investigated 
analytically by Mohazzabi [l(| • But no closed form for 
the quantisation is available in this case and the lack of 
regularity in the spacings has been computed numeri- 
cally. 

Finally, Stillinger and Stillinger p| have proposed a 
two-parameter family of isochronous potentials interpo- 
lating between the harmonic oscillator and the asymmet- 
ric parabolic well. Although not explicitly mentioned in 
their paper, this family includes the radial harmonic po- 
tential (also known as the isotonic potential) and Urabe's 
potential. For non extremal values of the parameters, 
this family is analytic on the entire real line. It also 
includes the one-parameter family of analytic potentials 
defined by Bolotin and MacKay [3] • To prove that the ex- 
act spectrum of these potentials is not equidistant in gen- 
eral, the authors based their argument on the assumption 
that energy levels depend continuously on the parame- 
ters. As the family interpolates continuously between the 
harmonic oscillator, whose spectrum is equispaced, and 
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the asymmetric parabolic well, whose spectrum is not eq- 
uispaced, they deduced that the spectrum is not exactly 
equidistant for some range of the parameters. However, 
this result is only qualitative in the sense that the au- 
thors do not provide a method to evaluate the spectrum 
for general values of the parameters. 

The aim of the present paper is precisely to fill in this 
gap and to provide a way to calculate the first corrections 
to the equispacing of the EBK spectrum. In sectionllTl we 
first review some properties of isochronous potentials and 
derive new analytical expressions for some one- or two- 
parameter families of isochronous potentials. We then 
use the perturbation method WKB up to fourth order 
(beyond the EBK quantisation) to derive an expression 
for the corrections to the equispacing valid for any ana- 
lytic isochronous potential (section lTTT|l . In section llVl we 
apply these results to the potentials described above and 
derive an analytical or asymptotic expression for their 
higher-order WKB quantisation condition. The resulting 
higher-order WKB spectra are then checked against nu- 
merical evaluations of the exact quantum problem and 
prove not only to reproduce the right behaviour of the 
quantum levels at high energy but also to be quite accu- 
rate at low energy for some range of the parameters. A 
discussion of these results and some conclusions are given 
in section IVll 



II. ISOCHRONOUS POTENTIALS 
A. Generalities 

The literature about isochronous potentials (or more 
generally about isochronous centres) is particularly vast 
and the interested reader is referred to |13j and refer- 
ences therein for more information. This section is not 
intended to present a detailed analysis of isochronism but 
merely to recall some basic facts and to introduce certain 
notations which will be useful in the rest of this pa per . 
Our approach of isochronism is mainly inspired by |l 1| 
and based on the so-called S function (see below for a 
definition). Another fruitful way to obtain explicit an- 
alytical expressions for isochronous potentials has been 
developed in [l5| . 

A potential V(x) is said to be isochronous if it gener- 
ates a motion x(t), obeying x + V '(x) = 0, whose period 
(frequency) is energy-independent. It can be shown that 
such a potential has a single minimum (see for instance 
0). For potentials which are at least twice differentiable 
at the origin, we will consider T^(0) = 0, V^'(0) = and 



V"(0) 



throughout, without loss of generality. 



For arbitrary potentials, the period of motion is 



x+(E) 



T{E) = V2 J 



dx 



c_(E) 



V E - V(x) 



(1) 



at energy E. Conversely, it has been shown by Landau 
and Lifshitz that, for a prescribed period T{E), there 
is a corresponding distance between the turning points 
given by 



x+(E) -X-(E) = 



T{u) 
rrV2j y/E~i 



du 



(2) 



Imposing isochronism to the potential (T(E) — est 
2-k/uj) yields 



x+(E) - x-{E) = 2- 



(3) 



If we let x(E) be the middle of the segment 
[x—(E),x+(E)], we obtain the equivalent formulation 



x±(E) = ±- 



>2E 



HE) 



(4) 



An entire freedom is left in the choice of x(E) provided 
the inversion of leads to a potential V(x) (= E) single- 
valued on the real axis. The smoothness of x induces 
the smoothness of V. Geometrically, the function x(E) 
defines a shear of the parabolic potential giving rise to 
another isochronous potential (for an illustration of such 
a transformation see 0). 

For our purpose, it's convenient to rephrase the results 
above in a slightly different way. Let us define a variable 
Igl through the transformation 



The two eqs. (QJ can now be recast into a single one 



(5) 



x{X) = X + x(X) 



(6) 



The function x(X) now represents the two branches of 
the potential V parametrised by X . The new shear func- 
tion x is related to x by x(X) = x(^uj 2 X 2 ) and is even 
in A, 



x(X) = x(-X) 



(7) 



As V(0) = 0, equations © and © yield x(0) = x(0) = 0. 
The condition dx/dX > of JSJ) ensures that X and x 
are in bijection, i.e. that V(x) is single- valued. Notice 
that we could have chosen dx/dX < instead. Now, 
from ©, 



where 



(8) 



(9) 



where x± (E) are the two turning points of the trajectory 



The latter is similar to the function introduced in [T^ 
but our definition of X does not include the frequency w. 
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Then, dx/dX > =>■ > -1. Moreover, as x(X) 

is even, S'(X) is odd, 



S{X) = 
which finally yields, 

VleR, |5(X)|<1. 



(10) 



(11) 



Remark that the function S defines the potential V up to 
a multiplicative constant only. Indeed, from JSJ we have 



x 



x = X + J S(u) du . 
o 

Inverting this relation gives X{x) and thus, from @ 
V(x) = ±u, 2 [X(x)} 2 . 



(12) 



(13) 



The multiplicative constant is thus basically the square 
of the frequency, which can be chosen arbitrarily. 



B. Scaling properties 

In what follows, we derive some simple, though 
important, results regarding the scaling properties of 
isochronous potentials. 

Claim 1. Let V(x) be isochronous with frequency to. 
Then 



V(7,/?) e 



V(x) = [j) V(J3x) (14) 



is isochronous with frequency Co = jlu. 

Proof. Let Iq(E) be the action associated to V(x). Then 

Jo (£) = ^0z]yjE-V{x)dx where V(x ± (E)) = 

E. Let Iq(E) be associated with V(x). Its turning points 
at energy E are related to those of V(x) by 



Hence, 



, . 1 f(3 2 E 
x±{E) = -x ± ^~ 



1 T fP 2 E 



I 0( E )=p2 I 0\ .o 



(15) 



(16) 



So far, this relation is true for all potentials. 

For an isochronous potential, however, Io(E) = E/u>. 
Then, from the previous relation, Io(E) = E/(^lu), which 
proves that V{x) is isochronous with frequency ju>. □ 



Corollary. // V{x) is isochronous with frequency u>, 

v(x) = jpV<M 

is isochronous with the same frequency. 



(17) 



This shows that, once an isochronous potential V(x) is 
known, a one-parameter family of isochronous potentials 
with the same frequency may be derived from the above 
scaling. 

Claim 2. The functions x(X) and S(X) corresponding 
to V(x) are related to x(X) and S(X) defining V(x) by 

x(X) = -=x{fiX) and S(X) = S(f3X) . (18) 
P 

Proof. By definition x(X) = X + x(X) which yields 
x{X) = \x{X) + x{-X)}/2. Starting from V(x) = 
\Co 2 X 2 , we obtain (7/ (3) 2 V{[3x) = ^ 2 lu 2 X 2 , that is 
V(/3x) = \uo 2 {(3X) 2 . Hence, x(X) = ^x(f3X). Rein- 
stating this last expression in x(X) and using x{X) = 
[x(X) + x{— X)]/2, we get the first of the relations 
iJTSJ. The second is obvious due to the relation S(X) = 
dx{X)/dX. □ 

Note that, as expected, (3 is the parameter involved in 
the scaling relations Ijl8|l while 7 is an overall prefactor 
responsible for the tuning of the frequency only. 

Claim 3. Asymptotic behaviour ofV(x) defined by j_?7| ) 
with respect to (3. Provided S(X) is defined on the entire 
real line, 



As 0^0, V{x)^- 



2 Vl + S +sgn(x) 



Asp^oo, V{x)^\( " 

2 \1 + (S)sgn(x) 



(19) 



(20) 



where 

S+ = lim S(X) and (S) = lim i / S(u) du . 

X^0+ T—*oc 1 Jq 

The function sgn(a;) is equal to 1 if x > and —1 if 
x < 0. If S(X) is continuous in X = 0, then S(0) = 0. 
The first expression then reduces to the harmonic oscilla- 
tor potential. The second one corresponds to a split har- 
monic oscillator whose frequencies are ui± = uj/(1 ± (S)). 

Proof. Let us rewrite equation ljl2|l in the form 

x = X ^1 + J 's{u)duj . 

Now, applied to the rescaled potential for which S(X) = 
S([3X) (see ©), this gives 



liX 



X \ l + W I s(v)dv 



Taking the limits (3 — > + and /3 — ► 00 yields x = 
X(l + S+sgn{X)) and x = X (1 + (S)sga(X)) respec- 
tively, and, as sgn(X) = sgn(ir), we obtain (|19fl and 

pn|i. □ 
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Now, using the results of the previous proof, we can 
obtain the asymptotic behaviour of an isochronous po- 
tential defined on the entire real line when \x\ — > oo. 

Claim 4. Providing \(S)\ ^ 1, where (S) is defined in 
claim\^ 



V{X) 2 (l + (S)sgn(x) 



x 2 + o(x 2 ), \x\ -> oo (21) 



When \ (S)\ = 1, the potential is singular. In this case, 
only one branch (say x — > oo) is asymptotically parabolic 

and such that V{x) = ^x 2 + o(x 2 ), x — > oo. 

C. Examples 

We now apply the above considerations to the determi- 
nation of isochronous potentials whose analytical expres- 
sion can be given explicitly. We start from the function 
S(X) which is required to satisfy the conditions given 
in ltTU|) and (JTTJ. We then apply scaling (|T7jl to derive 
a family of potentials with the same frequency. To our 
knowledge, the families of potentials derived in section 
III C 21 never before appeared in the literature. 



1. Family I 



Let 



S(X) 



aX 



Vl + aX 2 



a e 0, 11 



(22) 



The condition on a ensures that < 1. Moreover 

S(X) = -S(—X) so that, conditions JTOJl and (JTTJ) are 
fulfilled. The shear function corresponding to S(X) is 

x{X) = /* S(u)du = y/l + aX 2 - 1 and inverting O 
leads to 



X 



(x + 1) - ^ax(x + 2) + 1 



1 - a 

Using l|13[l we obtain the potential 



V(a;x) 



(x + 1) - y/ax{x + 2) + 1 



1 — a 



(23) 



This is precisely the one-parameter family of potentials 
derived by Bolotin and MacKay in Notice that the 
parameter a is distinct from the scaling parameter (3 in- 
troduced in the previous section. Thus, relation (|17f) en- 
ables us to derive a two-parameter family of isochronous 
potentials 



V(a,/3;x) = — ^ 



uP_ ( {fix + 1) - y/a/3x(0x + 2) + 1 



(24) 



presented for the first time by Stillinger et al. in pj (cf. 

The family (|24Jl includes several known isochronous po- 
tentials as limiting cases which are listed below. 



• a — ► or [3 — > 0: Harmonic oscillator. 

U 2 9 

V{x) = —x 2 . 

• a — > 1 and /3 7^ 0: Isotonic potential. 
,2 / , n 2 



(25) 



a; 



1 



/3a; + 1 



z > -- (26) 



This potential, also called radial harmonic oscilla- 
tor, has the property to be singular at x — — 1//3 
and to be defined on a half-line. 

(3 — ► 00: iSpMi harmonic oscillator. 

{if " "\ a;2 a; > 
2 ) 1+ ^( 2 ' " (27) 

In the case a — ► 1, the left part (a; < 0) of the 
potential converges to a hard wall in x = and 
gives rise to a half-parabolic potential. 



• a/3 = C 7^ and (3 — > 0: Urabe's potential. 

v2 1 3 



V(x) 



(28) 



Notice that in this last case, the condition a £ [0, 1] 
has been relaxed. Preserving |5*(X)| < 1 amounts 
to restricting the values of X to \X | < 1/|C|- Hence, 
a potential defined only on a finite interval of the 
real axis. 

Note that the first and third results can be obtained 
by means of the general relations (|19|) and 12011 . 



2. Family II 



Let 



g W * S ^ , (a, £) e M*+ 
£ cosh A — 1 + a 



(29) 



S'(X) is readily odd. According to the value of a, the 
condition IS^-X")! < 1 imposes 



< a < 1 => £ > [a(2-a)]- 1 /2 
a > 1 =*> £ > 1 



(30) 



The shear function of (|29|) is given by x(X) = 
{rn[(cosh(X) - 1 + a)/a]}/£ and inverting (T2J (or ©) 
amounts to solving 

+ 2(a- l)y« + F f ^ 1 - 2ae^ x = (31) 
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where Y = expX. For some particular values of a and 
£, the above algebraic equation can be solved exactly 
for Y. This leads to an analytical expression for the 
potential. In what follows, we work out some of the 
simplest cases and systematically apply the scaling (|17|) 
to derive the corresponding family. 



£ = 1 and a > 1. 



hi 



1 - a + ^2a{eP x - I) + a* 



x > 



-~ln2a. 



(32) 



This represents a two-parameter family of potentials sin- 
gular in x = — (ln2a)//3 and defined on a half- line. 
Notice that for a = 1, the analytical form of this po- 
tential simplifies further to yield the remarkably simple 
expression 




x > 



In 2 



1/3 



(33) 



(34) 



III. BEYOND EBK 

We now turn to the problem of the determination of 
the quantum spectrum of isochronous potentials. As al- 
ready stated in the introduction, their semiclassical EBK 
spectrum (i.e., obtained from a first order WKB method) 
is perfectly regularly spaced. Nevertheless, a study of the 
exact spectrum of the split-harmonic oscillator leads to 
the conclusion that at least not all isochronouspoten- 
tials possess strictly equidistant energy levels (0,0). As 
the determination of an exact quantisation condition is 
not possible in general, a natural idea is to go beyond 
the usual semiclassical approximation and to study the 
properties of the spectrum by means of the higher order 
terms generated by the WKB method. 

The WKB method to all orders has been first devel- 
oped by Dunham UM and subsequently improved by 
many authors 0, Il8l fill l2lj . In particular, it allows 
the quantisation condition for a ID analytic potential 
to be written as a power series in h. Such series have 
to be interpreted as asymptotic series and are generally 
not convergent. However, the result leads to the exact 
quantisation condition for those exactly solvable poten- 
tials, whose WKB series can be evaluated explicitly and 
summed. 

In what follows, we briefly recall the WKB method 
before to make use of the second and fourth terms de- 
rived in 01 to express the first corrections to the EBK 
quantisation in terms of the S function presented in the 
previous section. 



Contrary to the previous family, this one is defined on 
the whole line R and is not singular. 



WKB to all orders 



£ = 3 and a = 1. 



™ - w 



111 



VT+8e^- 1 



-1 2 



(35) 



is defined on the whole line R and is not singular, as in 
the previous case. 

Due to Cardano's formula, (|3TT) can be solved analyti- 
cally in the more general case where £ is either equal to 2 
or 3 and a > 1. For the same reason, it is also possible to 
derive exact solutions of l|31|) when a = 1 and £ = 5 or 7. 
These expressions, however, are somewhat messy and are 
not provided. It is also possible to derive exact solutions 
of H31JI when a = 2 and £ = 5 or 7. The potentials given 
above are a few among many others. 

Whatever the values of the parameters £ and a, the 
limiting case f3 — > leads to the harmonic potential as 
S(X) given by l)29|) is continuous in X = 0. Given that 
(S) = l/£, the limiting case (5 — > oo gives rise to a split- 
harmonic oscillator whose left and right frequencies are 
u)± = £u;/(£± 1). In the first case, where £ = 1, this 
leads to a half-parabolic well. 



For a complete description of the method and the 
properties of WKB series, the reader is referred to 
|17l llSl Y\S\ I2H . The present introduction is mainly in- 
spired by We start from the Schrodinger equation 
in ID (m = 1), 



2 dx 2 



+ V(x) 



tp(x) = Eip(x) 



Writing 



we obtain 



ip(x) = exp ( ) , 



a' 2 (x)+^-)a"(x) = 2(E-V(x)) 

solved by means of a power series in h, 

k 



*>-E(?) 



a k (x) . 



(36) 



(37) 



(38) 



(39) 
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Reinstating (|39() in (|38|l and solving order by order in h 
yields the recurrence 

I 

a' 2 = 2(E-V) and jV^oU + ^-x = 0, / > 1 (40) 

Requiring the wave function to be single valued amounts 
to imposing the following quantisation condition 



dak = 2nhn , n G N 



(41) 



where the complex integration contour 7 surrounds the 
two turning points of V(x) at energy E located on the 
real axis. 

The first term of the series H41JI is readily proportional 
to the classical action Iq(E) 

-(E) 



da a = 2 / y/2(E-V(x))dx = 2tt I (E) . (42) 

As integral of a logarithmic derivative, the second term 
can be shown to be equal to 



dai = —Trh . 



(43) 



If we truncate the WKB series to this last order, we re- 
derive the so-called EBK quantisation 



Io(E) = 



1 



h, n G 



(44) 



where the 1/2 term represents the Maslov index. 

As shown by Froman |23j, all odd terms 02/s+i> ^ — ^ 
are total derivatives and as such, their contribution van- 
ishes. This allows us to rewrite the quantisation condi- 
tion as 



fc=0 

where we have defined 

hk(E) = 1- 
Zir 



21, 



h, n G N 



d<r 2k , k G N . 



(45) 



(46) 



In case V(x) is analytic and V'(x) ^ if x ^ 0, the 
authors of [21j have proved that the contour integrals 
involved in (|46|l can be systematically replaced by equiv- 
alent Ricmann integrals between the two turning points. 
Thanks to this result it can be shown that 0] 



h{E) = - 



x+(E) 

d 2 f , v 2 {x) 



dx ■ 



24V2irdE 2 J JE~^ V (x) 

X -{E) 



(47) 



and 



h(E) 



4V2tt 



-(E) 



dx ■ 



V" 2 {x) 



1 c> 3 

UOdEP J ~ JE - V(x) 

X -(E) 



x+(E) 

id 4 r , v 2 (x)v"{x) 



dx ■ 



288 8E* J ^JE - V(x) 

x.(E) 



(48) 



As the H's are originally given by contour integrals, 
addition of total derivatives to the a' 2k doesn't change 
their value. However, it modifies their expression in 
terms of the potential and its derivatives and possibly 
allows for their simplification. Such a technique is devel- 
oped and used systematically in [21[. Expressions <|47Jl 
and H48fl are one among others (see for instance po|). 



B. WKB expansion for isochronous potentials 

We now use the results obtained in section III Al to ex- 
press I2 and Z4 in terms of the function S(X) related 
to the potential V(x). We make use of the change of 
variables V(x) = uj 2 X 2 /2 and of the relation JHJ) to find 



dV uj 2 X 
~dx ~ l + S(X) 



(49) 



Remarking further that X(x±(E)) = y/2E/u>, we obtain 



h 2 ^ d 2 



X 2 dX 



24V2irdEi J ^ + S (X))Je - \u*X* 



h{E) 



(50) 

Then using the fact that S(X) is odd and making the 
change of variables u = u 2 X 2 /2, we find 



HE) 



h 2 m d 2 

12^'dE 1 



E 



,1/2 



du 



(1 - u) 1 / 2 1 _ g2 ^VjEhL^j 



(51) 

Another equivalent formulation suitable for an asymp- 
totic analysis is given by 



HE) 



h 2 u> 1 

12^ & J {E-vy/ 2 dv 2 




v 



dv . 



This way, I 2 is expressed through an Abel type (or 
Riemann-Liouville fractional) integral about which a lot 
is known [2II3I. 
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A similar calculation yields the fourth order correction 



h(E) 



1 r< ( y/2Eu \ 

1 d 3 f G i v^^J 



120 9£; 3 J U 1 /2(l_ u )l/2 




du 



t V2E 7 



i d 4 

288 ae 4 7 M 1 /2(i_ u )i/2 





(53) 



where the functions G\ and G2 are given by 

• 4 r 8JC SS'(1 + S 2 ) 



Gx{X) = 



G 2 (X) = 



3S^ + 1 



(1 - S 2 ) 3 
A 2 S' 2 (l + 10S 2 + 5S 4 )i 



1-S 2 



(1 - S 2 ) 2 



LO b X 2 



(1 - S 2 f 



3S 2 + 1 



AXSS'{\ + S 2 ) 
1-S 2 



(54) 



In the last expressions, S stands for S(X) and 5" for 
dS/dX. 

Similar to 1% , I4 is equivalently expressed through Abel 
integrals via 



HE) = 



A-kuo 



E- 



E -3 J y 5/2 d 3 

120 J (E - v) 1 / 2 !^ 





V 7 ' 2 



288 J (E-v^dv 4 




G 2 




C. Scaling properties of WKB series 



We have seen in section III Bl that the general scaling 
V(x) = (-f/f3) 2 V(/3x), (7,/?) G R* 2 , preserves isochro- 
nism and changes the frequency u> of V(x) to u> = 70;. 
Such a transformation, however, is not expected to have 
simple consequences at the quantum level given the lack 
of relation between the spectra of V(x) and V(x) in 
the general case. This can be seen directly from the 
Schrodinger equation itself 



-Ydx^ + V{x) 



2 dy 2 



Ei{j(x) (56) 
E 

w 



Mv) ( 57 ) 



where y = f3x and ip(y) = ip( x )- 

This last equation is equivalent to the Schrodinger equa- 
tion for V iff 7 = ±/3 2 , which corresponds to the particu- 
lar scaling V{x) — /3 2 V(/3x). The spectrum of V is then 
the spectrum of V multiplied by 1 . 

But, nothing can be said about the general trans- 
formation and the frequency- preserving scaling V(x) = 
V(f3x) I p 2 in particular without further information re- 
garding the potential V(x). 



As we are dealing with isochronous potentials, we 
may transform the Schrodinger equation 13(j|) into an 
equivalent Sturm-Liouville problem involving the func- 
tion S(X) rather than the potential V(x) itself. Thanks 
to JSJ) and ||SJ|, we obtain 



Ydx 



1 



1 + SdX 



+ (1 + 5) 



E 1' 

2 



= (58) 



where 4> = <j>{X) = ip(x) and S = S(X). This last form 
proves convenient for numerical purposes as S(X) is a 
supplied function whose analytical form can be chosen 
freely within the requirements (|10|l and (|11|) . But we 
found it of no particular help in investigating scaling ef- 
fects on the spectrum. 

Instead, we can look at the way this scaling affects each 
term of the WKB series. 

Claim 5. Letl-x n andl 2n he the 2n th terms of the WKB 
series for the potentials V(x) and V(x) — (7 / '/3) 2 V \(3x) 
respectively. Then, 



hn(E) 



2n-l 



'2)1 



(3 2 E 



n e N. 



(59) 



Notice that this expression is valid for any potential 
and not only for isochronous potentials. Details of the 
proof are given in appendix^ But let us immediately re- 
mark that, multiplying equation (|57|l by (3 V7 2 indicates 
that the energy is rescaled by (3 2 /^ 2 and h, by /3 2 /7- 
Since each term of the WKB series is proportional to a 
given power of h, the above result is easily established. 

Thus, we arrive to the conclusion that, even though no 
simple transform enables us to derive the spectrum of V 
from the spectrum of V, each term of the WKB series 
obeys its own scaling, which allows for calculating the 
WKB series for V once it is known for V. We may restrict 
our investigation to potentials with a fixed frequency and 
a fixed value for [3. 



IV. APPLICATION 

We now apply what precedes to some isochronous po- 
tentials presented in section |nj We start our investiga- 
tion with the family of potentials derived by Bolotin and 
MacKay and obtain an analytical expression for its 
fourth order WKB quantisation condition. This result is 
immediately generalised to the two-parameter Stillinger's 
family (family /) thanks to the scaling (|59|l . 

Regarding the class of potentials presented in III C 21 
(family II), we won't be able to derive any explicit ana- 
lytical expression for the corrections. We will neverthe- 
less derive their asymptotic behaviour. 



A. Family / 

For the sake of simplicity, we calculate the terms I2 and 
I4 for a potential of frequency u> — y/2 whose S function 
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is given by 



S(x) = v/T^ 



X 



Vi + x 2 



, r,e [0,1]. (60) 



We recover a general expression of the type l|24l) for the 
potential by means of ifTHl) with S(X) = S(f3y/aX) and 
77 = 1 — a. The frequency is then fixed to the value ui 
thanks to 7 = u/ \/2. 



1. Analytic expressions for I2 and 1 4 

Setting h=l, expression l(5*T|l leads to 
V2 (1 - r?) 



h{E) = - 



2 4 (I + 77.E) 5 / 2 



and expression (|53"j) to 



V2 (l-ry)pfag) 
2io (1+r/S) 15 / 2 



(61) 



(62) 



where 
p{rj, E) 



l5r] 4 E 4 - 30ri 3 {7r) - 16)E 3 
3f7 2 (119?7 2 - 651?? + 411)£ 2 
3r?(455?7 2 - 567?7 + 104) E - 280? ? 2 + 140?? + 8 



2. WKB quantisation condition to fourth order 



in this limit. It can be shown by calculating the entire 
WKB series that all corrections are energy-independent 
and that their summation leads to the exact quantisation 
condition 0, [2^ • We can check that the two corrective 
terms I2 and I4 we have obtained are the coefficients of 
the Taylor expansion for the non-integral Maslov index 
(see |24j) of the exact quantisation condition in j3 = 0. 
Indeed, the exact quantisation condition for the radial 
harmonic oscillator i|26[l (or isotonic potential) reads (see 
for example [25|) 



E n 



(n + -)u 



where the non-integral Maslov index is given by 




(65) 



Expanding this last expression around (3 = yields 



1 1/3 2 1 f3 6 

E n ~ (n+ - + — -r 

2 8 w 32 w 3 



which is precisely what we obtain from (|63|l given that, 
a = 1 implies Q = 1 and P = 8. 

This remark brings us to consider that WKB series 
should be more and more accurate as the scaling param- 
eter (3 tends to zero, i.e. in the limit where the potential 
tends to the harmonic oscillator. This statement is con- 
firmed by the numerical computations performed in the 
next sections. 



The fourth order WKB quantisation condition (|45|l for 
the family of potentials (I24|) is deduced from (|61|) and 
(|62|l and from the scaling relation l|59l) 



E- 



1 (a/3) 2 a 4 /3 6 P 



8 Q5/2 2 8 w 2 Q 1 5/ 2 



(63) 



where 



1 



I3 2 ' 

2a(l-a)-^E 

UJ 2 



a 2 



P = p [l - a,2a—E 

(64) 

We can now solve equation H63J1 for E to obtain a 
semiclassical approximation of the spectrum. This ap- 
proximation should be valid asymptotically, that is, as 
E — ► 00, as we expect I (E) < h(E) < HE)---. 
We can check that this is indeed the case for a < 1 
as HE) cx E, I 2 {E) oc E^' 2 and h(E) cx E- 7 ' 2 . No- 
tice that for the particular value a = (or (3 = 0), the 
potential l|24|l becomes a harmonic potential and the cor- 
rections I2 and I4 vanish. It can be shown that all higher 
corrections Ii n vanish as well 0|i which establishes the 
well known fact that, in this case, the EBK quantisation 
leads to an exact result. 

However, for a = 1, the potential l|24|) becomes a sin- 
gular isotonic potential (see section III C 1|) . The sec- 
ond and fourth corrections become energy-independent 



B. Family II 

We restrict our study to the nonsingular potentials l|34() 
and l|35|) for which the WKB method leads to an asymp- 
totic series as E — * 00 40]. In what follows, we again 
evaluate 1% and I4 for lu = y2. The S function of lp?4^ 
and reads 

St(X) = j tanh(X) , Ce{2,3}. (66) 

We recover the general expression for the potentials by 
means of the scalings (|14|) and l|18fl . The frequency is 
brought back to the general value u> by choosing 7 = 



1. Expressions for I2 and 1 4 

Although S^(X) in l|66|l has a very simple form, we are 
not able to provide an analytical expression for I2 (E) and 
14(E). Then, we resort to their numerical evaluation as 
a function of E and use the corresponding functions to 
calculate the semiclassical spectrum. 

The lack of exact analytical forms for I2 (E) and I4 (E) 
does not prevent us from extracting some useful infor- 
mation regarding their asymptotic behaviour as E — > 00. 
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TABLE I: Table of the coefficients M m ,n(£) appearing in 16 7H 
for uj = y/2. 





£ = 2 


C = 3 


M 2 ,i 


f .9020 x 10~ 2 


6.3076 x 10~ 3 


M 2 ,2 


f.7287 x 10 -1 


5.5608 x 10~ 2 


M 4 ,i 


7.0f28 x 10~ 3 


f.94f4 x 10~ 3 



lutely essential. It comes from the fact that, when deriv- 
ing the asymptotic forms i|67[) . we have assumed that E 
was large in h(E) and Ia{E). As fflfy rescales the energy 
by a factor (f3/uj) 2 , we have to consider that {J3jui) 2 E is 
now large. For example, this condition ensures that the 
corrective terms of order E~ 5 / 2 and E~ 7 1 2 tend to zero 
in l|68|) as (3 — > (harmonic limit) providing E > (3~ 2 . 



According to the fact that WKB should provide the best 
results in this limit, we expect this analysis to yield the 
right correction to the harmonic levels at high energy. 
Using H52f> and 155fl and a result regarding the asymp- 
totic expansion of fractional integrals obtained in 27] (see 
also |41|), we easily derive the large-energy behaviour of 
h{E) and h{E) to be 



M 2 ,i(Q , M 2 , 2 (Q 

£5/2 

M 4 ,i(Q 
EV 2 



+ o 



E 7/2 

1 



E 7/2 



(67) 



where £ £ {2,3} and M mj „(£) are some nonzero coef- 
ficients related to the Mcllin transform of the functions 
appearing in the numerator of the fractional integrals in- 
volved in l(5^|l and JHSJj (see or [2]} for more details). 
The numerical values of these coefficients for a frequency 
uj = \/2 are listed in table [I] 



2. Asymptotic quantisation condition 

Our asymptotic expansion in E is truncated at order 
E~' 1 1 2 in l|67|l because the term Iq{E) would be needed 
to evaluate the proper coefficient of order E~ 9 / 2 for the 
entire WKB series. We note again on this example that 
WKB series are generally no direct asymptotic expan- 
sions in energy. Obtaining such a series requires l2 n (E) 
to be expanded as asymptotic series before terms of equal 
power in E can be collected. A method has been devel- 
oped in [2^| which allows for the derivation of a proper 
(direct) asymptotic expansion in E. Its range of applica- 
tion seems to be restricted to power-law potentials, how- 
ever, and it is not clear how it could be used in the present 
case (see |29ll30|l. 

Finally, applying the scaling to i|(77|) . we obtain 
an asymptotic quantisation condition for potentials i|34|) 
and l|35|) valid for any frequency u> and any parameter j3 
such that (f3/uj) 2 E is large 



E 



uj 5 M 2 ,i 

4[3 3 E 5 / 2 ' 4/3 1 2/3 



-713 i 7T3I M 2,2 + M 4 ,i 



£7/2 



(n+ 2 )w ' ^ 



(68) 



For simplicity, we have dropped the £ dependence of the 
coefficients M mjn (£). The last condition in l|68|) is abso- 



V. NUMERICAL EVALUATION OF THE 
SPECTRUM 

A. Numerical solution of singular Sturm-Liouville 
problems 

As the difference between the spectrum of a general 
isochronous potential and the harmonic one is expected 
to be small, its numerical determination requires an ac- 
curate method. 

Our numerical evaluation of the spectrum is based on 
the Sturm-Liouville problem - SLP - 158|) . We have solved 
this equation by means of a shooting method provided 
by the code SLEIGN2 developed by Bailey, Everitt and 
Zettl [32Ll33T |. Among other advantages, solving l|58|l only 
requires the analytical expression for S(X) and not the 
potential. It can then be easily used for a wide range of 
potentials. 

The SLP (|58|l is always singular given that its end- 
points are located at X — ±00 (where the function (j)(X) 
vanishes). The type of singularity may vary according to 
the potential under consideration. But in any case, it is 
the same for the SLP H58fl and the associated Schrodinger 
equation l(S^|) (see [33 lemma 1 and 2). 



B. Endpoint classification 

Runing SLEIGN2 requires knowledge of the singularity 
type of each endpoint of the SLP under investigation. 
We, then, classify hereafter those susceptible to be met in 
handling isochronous potentials. To determine the type 
of singularity of each branch of V(x), we first remark that 

Claim 6. IfV(x) is isochronous of frequency to 



Vx e 



V (x) > ^-x 2 



(69) 



This is a consequence of the fact that (S'(A')I < 1 
\x(X)\ < \X\. Now, x = X + x{X) =*> |x| < 2\X\ hence 
u?x 2 l% < uj 2 X 2 /2 = V(x). 

Claim 7. An endpoint x e of V(x) located at ±00 is al- 
ways of both the limit-point and nonoscillatory type. 

Indeed, applying theorem 7 of |35| and claim [5] proves 
this endpoint to be of the limit-point type. Moreover, 
as \mi x ^ Xii V{x) = 00, applying theorem 6 of the same 
paper shows that x e is of the nonoscillatory type. 



10 



Claim 8. Assume an endpoint x e of V(x) is located at 
Xq, where Xq is finite. Then it is always of the nonoscil- 
latory type. 

• // lim inf (x - x ) 2 V(x) > 3/8 or 

X — >Xq 

if(x-x ) 2 V(x) > 3/8 and lim inf (x - x Q ) 2 V(x) = 3/8 

X — >Xq 

then, the singularity is of the limit-point type. 

• If lim sup (x - x ) 2 V(x) < 3/8, 

x — >xo 

it is of the limit-circle type. 

Indeed, as liirLc—^o V(x) = +oo, applying theorem 4 of 
|35j shows immediately that the type of singularity of the 
branch located at xq is always nonoscillatory. Whether it 
is of the limit-point or limit-circle type relies on theorem 
5 of the same paper. 



C. Family I and Family II potentials 

Now, we investigate the spectrum of real-analytic po- 
tentials belonging to Family I or II numerically and com- 
pare it to the WKB predictions. We select two specific 
types of potentials for this study with a frequency set to 
unity. 

For potentials of type I, we consider Ij24(l with a = 1/2 
(a < 1 implies that the potential is nonsingular) and we 
denote it 



K3 



l 

0.9 
0.8 
0.7 



H|n 0.6 
+ 



I 



0.5 - 

0.4 - 

0.3 - 

0.2 - 

0.1 - 




Exact + 
WKB 







30 

E,, 



FIG. 1: Comparison of the exact and semiclassical corrections 
to the harmonic spectrum for the first 50 energy levels of 
V/(l/2; a;) (see (I7UH '). Energy differences e n and e WKB (E) (see 
17211 and 1741 1 ) are multiplied by E 5 ^ 2 . Cross symbols are the 
numerical results and the solid line indicates the WKB results. 



We define the exact (numerical) correction to the har- 
monic levels by 



e n = E n - (n + 



(72) 



Vi(P;x) 



f3x 



((3x + 2) + 1 



and, according to (|45|l and l|59|l . its (continuous) fourth 

(70) 

order semiclassical approximation by 



Of course, some quantitative differences are to be ex- 
pected according to the value of a but, qualitatively, 
the features we are going to discuss hereafter remain un- 
changed. 

For potentials of type II, we select ().'.'> ■"> I) and denote it 



V n (i3;x) 



1 



In 



Vl + 8e 3 ' 3x - 1 



(71) 



Again, we could have chosen (|34|l instead without conse- 
quences regarding the qualitative behaviour of the spec- 
trum as a function of the parameter (3. 



1. Small f3 regime 

As already noted on the specific example of the iso- 
tonic potential (see section HV A 2|) . the WKB method is 
expected to provide the best results in the limit where the 
parameter (3 is small. In this limit, the potential is close 
to being harmonic. We verify this statement by compar- 
ing the exact numerical corrections to the harmonic lev- 
els to those obtained from the semiclassical quantisation 
condition. For Vi((3; x), we use the analytical expression 
(|63|l whereas for Vjj((3; x), we resort to a numerical eval- 
uation of 12(E) and Ia(E). Within this section, [3 = 1/2. 



-0 1 h(f3 2 E) 



l3 6 h(P 2 E). 



(73) 



In this last expression, 12(E) and 14(E) are given by l|52() 
and (|55l> , respectively, and are evaluated at cu — 1 for the 
functions S(X) defined in {23 (type I) and (type 
II). For Vi(/3;x), §^ yields 



XE) 



1 S3 2 (3 6 



P 



2 5 Q5/2 2 12 Q15/2 



(74) 



1/2 



where P and Q, defined in (|64(l . are taken at a 
and uj = 1. 

The results are presented in Figures and For the 
sake of visibility, in both cases, the correction e is mul- 
tiplied by the inverse of its asymptotic decay, E b / 2 (see 
(jrjlfl and JHSJO- Cross symbols represent the exact nu- 
merical results and the solid line, the fourth-order WKB 
results. As we can see, the fit is virtually perfect for the 
two spectra. In the limit where f3 is small, WKB provides 
corrections which are very accurate even in the lowest 
part of the spectrum. And this is so even though, in 
both cases, E 5 / 2 e mKB (E) is far from reaching its asymp- 
totic values given by V2 (type I) and -2M 2 ,i ~ -1.261 
(type II), respectively. On these precise examples, we see 
that WKB carry much more information than a simple 
high-energy asymptotic expansion. 
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FIG. 2: Comparison of the exact and semiclassical corrections 
to the harmonic spectrum for the first 50 energy levels of 
V n (l/2;x) (see <|7TJ) multiplied by E 5/2 as in Fig. □ Cross 
symbols are the numerical results and the solid line, the result 
of the numerical evaluation of 1731 . 



FIG. 4: Comparison of the exact and split-harmonic correc- 
tions to the harmonic spectrum for the first 50 energy levels 
of Vn(30;x) (see 171l l). Like in Fig. |3] corrections are mul- 
tiplied by E 2 . Cross symbols are the numerical results and 
circles, the corresponding split-harmonic correction. 
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FIG. 3: Comparison of the exact and split-harmonic correc- 
tions to the harmonic spectrum for the first 50 energy levels 
of Vi(50;x) (see l)7()|l 'l. Notice that, contrary to Fig. Q cor- 
rections are now multiplied by E 2 . Cross symbols are the nu- 
merical results, circles are the corresponding split-harmonic 
correction and the dotted line is the analytical approximation 
Xn obtained in 1B3II (last expression). 



2. Large [3 regime 

As shown in claim in the limit of large f3, an 
isochronous potential asymptotically converges towards 
a split-harmonic oscillator. Thus, we expect the cor- 
responding quantum spectra to be close. For this rea- 
son, an analysis of the spectrum of the split-harmonic 
potential has been done in appendix [5] where a high- 
energy asymptotic expression for its quantum levels is 
derived. Using these results, we now compare the spec- 
trum of potentials Vj(50;x) and Vn{30;x) to the spec- 
trum of their respective asymptotic split-harmonic po- 
tentials (/? — > oo). 



According to l|27|) . when (3 — ► oo, Vi((3; x) converges to 
a split-harmonic oscillator with left and right frequencies 
whose ratio is given by p r — (1 — yfa)/(l + ^/a) = 3— 2y2- 
Regarding Vn(f3;x), this ratio becomes p H =1/2 (see 
end of section llV Bfl . 

Again, for the sake of clarity, we have multiplied l|72|l 
by the inverse of its split-harmonic asymptotic decay. 
Notice that the latter is different from the WKB decay. 
Indeed, as shown in ()B2f) . the split-harmonic correction 
to the harmonic levels vanishes like E~ 2 as E increases 
whereas the WKB correction decays like E~ 5 / 2 instead. 

In FigureEl even though e„ behaves erratically, in con- 
trast to the regularity observed in Figure^ the exact nu- 
merical results for Vi(50;x) are very well approximated 
by the split-harmonic corrections for low-energy levels. 
The apparent complex behaviour of e n is easily under- 
stood from expressions l|B2|) and l|B3(l of appendix [B] 
which show that, once multiplied by E 2 , the correction 
(Xn) i s a periodic function of E n (or n + 1/2) sampled at 
an incommensurate frequency. The dotted line of Figure 
13 which is a continuous version of \n (last expression of 
(|B3|) ). shows how accurate this asymptotic result is. 

In Figure 01 we have reported the results obtained for 
Vn (30; x). For such a large value of (3, the correction 
to the harmonic levels is again very close to the split- 
harmonic correction. But this time, due to the rational 
value of p H = 1/2, (|B3|) shows that E 2 e n ~ x« takes 
on two values only, ^toj and — in perfect agreement 
with those observed in Figure g| 

Yet, both Figure[3]and Figure0]indicate that the exact 
correction e„ and its asymptotic approximation progres- 
sively split up as the energy increases. This is what is 
expected if the WKB analysis becomes exact as E — > oo. 
Then, a crossover should exist between the two different 
power-law decays, E~ 2 , which is valid at low energy and 
i?~ 5 / 2 , which is valid at high energy. This is difficult to 
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ergy which are clearly reminiscent of those produced by a 
split-harmonic oscillator. It can be checked that they are 
almost of the same frequency as those observed in Figures 
Eland 01 But at higher energy, they are damped out and 
converge towards the asymptotic behaviour predicted by 
the semiclassical approach. Indeed, from (|7^jl and ll68l) . 
we easily obtain E 5 ^ 2 e VJKB — > V2/64 for Vi(2;x) and 
£ 5 / 2 e WKB -> -M 2 ,i/8 for Vn(^2;x). Hence, the asymp- 
totic values 2.21 x 1(T 2 and -7.88 x 1(T 4 oberved in 
Figures and respectively. 



VI. DISCUSSION AND CONCLUSION 



FIG. 5: Comparison of the exact and semiclassical corrections 
to the harmonic spectrum for the first 100 energy levels of 
Vi(2;x) (see I7UH 1. Energy differences e n and £ WKB (E) (see 
1721 and 1741 1 are multiplied by E 5 ^ 2 . Cross symbols linked 
by a dashed line are the numerical results and the solid line 
indicates the WKB results. 
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FIG. 6: Comparison of the exact and semiclassical corrections 
to the harmonic spectrum for the first 100 energy levels of 
V n (^2;x) (see J7TJ0 multiplied by E 5/2 as in Fig. Cross 
symbols linked by a dashed line are the numerical results and 
the solid line, the result of the numerical evaluation of l|73|l . 



observe for large values of (3, however, because it would 
require a very high numerical accuracy when going up 
the spectrum. Instead, we try to observe this transition 
for an intermediate value of the parameter /3 in the next 
section. 



3. Intermediate f3 regime 

In this last section, we select the two potentials V/(2; x) 
and Vn ( v^; x) . The results are shown in Figure |5] and 
Figure respectively. The correction e n is again multi- 

5 /2 

plied by the inverse of its asymptotic WKB decay, E n ' . 
In both cases, we observe some oscillations at low en- 



In this paper, we have provided a quantitative way to 
analyse the spectrum of analytic isochronous potentials. 
As is already known, the spectrum of an isochronous po- 
tential is generally not strictly regularly spaced, in con- 
trast to the harmonic one. Therefore, a method is re- 
quired to obtain the possible corrections to the harmonic 
equispacing of its energy levels. 

We have shown that, a WKB analysis going beyond 
the semiclassical EBK quantisation is precise enough to 
successfully account for these small differences. It pro- 
vides the right asymptotic behaviour of these corrections 
at high energy, and in certain cases, is even able to de- 
scribe them accurately from the ground-state. 

Moreover, we have seen that, for any isochronous po- 
tential V(x), there exists a simple scaling transforma- 
tion, namely V{f)\x) = V(f3x)/ (3 2 , which preserves its 
isochronism as well as its frequency. Under certain reg- 
ularity conditions, the one-parameter family V(/3; x) in- 
terpolates continuously between the harmonic oscillator 
(/3 — * 0), and a split-harmonic oscillator (/3 — •> oo). Gen- 
eral spectral features, with respect to the parameter /3, 
have been sketched on several examples. 

When is small, that is, when V(f3; x) is close to a 
harmonic oscillator, a higher-order quantisation condi- 
tion derived from the WKB series accurately describes 
the entire spectrum, even its lowest part. When (3 is 
large, V(f3; x) is close to a split-harmonic oscillator and 
the lowest part of its spectrum essentially reproduces the 
split-harmonic features. Nevertheless, at high energy, the 
spectrum becomes asymptotically described by a higher- 
order WKB quantisation condition again, which signifi- 
cantly differs from the exact split-harmonic spectrum. 

Corrections to the harmonic spectrum, that is, e n = 
E n — (n+^)huj, have been shown to scale like E~ 2 at high 
energy for a split-harmonic oscillator. Those generated 
by a higher-order WKB method, however, scale gener- 

— 5/2 

ically like E n . Assuming that the latter ultimately 
represent the right high-energy asymptotic behaviour of 
the spectrum, this leads to the conclusion that a crossover 
between these two trends exists. Although its precise lo- 
cation within the spectrum is difficult to explore numer- 
ically, the spectrum of V(/3; x) for intermediate values of 
the parameter /3 supports evidence of such a transition. 
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Notice finally that, in the examples treated in the last 
section, the first term of the WKB series going beyond 
the EBK quantisation, namely I2, dictates the leading 
order of the asymptotic decay of the correction e„ as 
E n — > 00. Moreover, in these examples, e n ~ E n 5 ^ 2 . 
We would like to point out that this needs not always be 
the case. Indeed, starting from expression (|52|) for I?{E ) 
and using the properties of Abel transforms |26| . |36j|. 
its possible to invert the problem and to calculate the 
general expression of the function S(X) corresponding 
to a prescribed function 12(E) (see appendix Q. 

Using this result, we can choose 12(E) (and deduce 
the corresponding analytic isochronous potential via the 
function S(X)) such that, its asymptotic decay is faster 
than the asymptotic decay of I±(E). An example given 
in appendix El proves that it is possible to find analytic 
potentials such that I 2 (E) ~ E~ 9 / 2 and h(E) ~ iT 7 / 2 
as E — > 00. The asymptotic leading order is then given 
by I4, in this case and different from the generic E~ 5 / 2 
decay. This precise example indicates how difficult it is to 
draw some general conclusions regarding the behaviour 
of the WKB series with respect to the energy. 

Our last remark concerns the class of isochronous po- 
tentials with a strictly equispaced (harmonic) spectrum. 
We have met such a one-parameter family in section 
IIV A 21 called isotonic or radial-harmonic oscillator. The 
interesting question is whether this family is the only one 
to be both classical and quantum "harmonic" . Without 
this being a proof, we would like to mention the following 
interesting result: 

Claim 9. The most general family of analytic 
isochronous potentials, such that all the terms of the 
WKB series, l2n(E), n > 1, defined in (|46|) . are con- 
stant (energy-independent), is the family of the isotonic 
oscillator with potential 

Proof. As already proved in [T^. |22| . all the terms 
l2n (E) , 71 > 1 of the isotonic potential are constant and 
once summed, the WKB series leads to the exact quan- 
tisation condition for this potential. To prove that this 
family is the most general, we use the result of appendix 
10 which shows that requiring 12(E) to be constant im- 
plies that S(X) corresponds to the function of the iso- 
tonic oscillator. □ 

Obviously, the fact that the I 2 n's, n>l, are constant 
ensures that the spectrum is strictly equispaced. Indeed, 
knowing that Iq(E) is proportional to E, because the 
potential is isochronous and that all other terms (I\ in- 
cluded) are constant, immediately leads to a quantisation 
condition of the form E = a(n + /i), (a, /j) 6 R 2 , and 
hence, to a regular spacing between consecutive energy 
levels. 

Unfortunately, the claim above is no proof that the 
isotonic family is the only one for which all terms of the 



WKB series add up to an energy-independent expression. 
For example, we could think of an overall cancelation of 
the energy-dependent terms generated by the J2 Tl (-E)'s. 
Moreover, although the last condition is quite a natural 
way to ensure that the spectrum is strictly equispaced, 
the function Y^=i^2n(E) could be energy- dependent 
and the quantisation condition X)^°=o I?n(E) = (n + ^)h 
still have an equispaced spectrum for solution. Although 
they seem quite unlikely to us, these possibilities cannot 
be excluded and the question of the class of potentials 
both classical and quantum harmonic remains open. 
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APPENDIX A: SCALING OF THE WKB TERMS 

Here, we show the transformation affecting each term 
of the WKB series as the potential V(x) is rescaled to 
V(x) — (7 / (3) 2 V((3x). First we consider the Schrodinger 
equation for V(x), 

h 2 d 2 ib(r) ~ 

"Y ~W~ + = E ^ • ( A1 ) 

According to (|45|l . its WKB quantisation condition reads 



J2lME)=[n+~] h, neN (A2) 
fc=0 ^ ' 

Given equation 1)46(1 . l2k(E) is proportional to h 2k and 
we will explicitely write it 

hk(E) = h 2 % k (E) . (A3) 

Let y — (3x, ip(y) = ip(x). Eq. (|A1|) is immediately 
transformed to 

4(^) 2d ^ + v(yWy) = my) (A4) 

where 

E=[^jE. (A5) 

Equation l|A4J> is a Schrodinger equation in the poten- 
tial V(x) with an effective Planck's constant h c g — 
Therefore, its quantisation condition is 

g(^) 3 W(„+i)^ ,„ eN (M) 
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Given that i|Al|l and i|A4|l are one and the same equation 
their quantisation conditions are the same. Dividing (IA6I) 
by — and identifying the series term by term to i|A2[l . 
we obtain the desired scaling 



hk{E) = ( 



02 \ / ;2 



I 2k 



(5 2 E 
7 2 



(A7) 



Another way to proceed is to consider the explicit ex- 
pression for l2n(E) given in [2l|. According to Eq. (44) 
of this reference, for m > 1, l2m(E) can be written as 

m 

h m {E)--— ^ (2m _3 + 2H)!! (A8) 



where 



Mm = -L^ : T I ^ (A 9) 



and where ^ = (y\, v%, . . . , i^m); € Z + , 
Y,]=i3 v j and |H = Si=i ^j- Moreover, 



(A10) 



Coefficients E/„ are defined by a reccurrcncc relation but 
are not given here since they are not needed in deriving 
the result. 

Let us evaluate l2 m (E) which is related to the potential 
V(x) = (-//(3) 2 V((3x). We first see that 

/ \2M 

^ W) ^ = \j) P L{v)v{v) {M- (All) 

Taking this result and l|15fl into account, we obtain 

J M (^) = i(^ 2(m " 1) ^)j M (^) (A12) 

Given that in Eq. I|A8J) . the sum is restricted to terms 
such that L(v) — 2m, we finally reobtain (|A7|) . 



APPENDIX B: SPECTRUM OF THE 
SPLIT-HARMONIC POTENTIAL 

The split-harmonic potential is made of two half- 
parabolic arches connected in x = with different fre- 
quencies u>i and uj r to the left and to the right respec- 
tively. It is an isochronous potential with frequency 
u = 2ujiuj r / '(u>i + uj r ). It turns out that an exact ex- 
pression for its quantisation condition is available and 
reads (see for example 0) 



1 



= (Bl) 



In this expression p = u>i/uj r is the ratio between the 
left and right frequencies and without loss of generality, 
we restrict its values to the range [0, 1] (values within the 
range [1, oo] give the same spectrum because they amount 
to flipping the potential around its vertical axis). Finally, 
x = (y + l/2)/(l + p) and E v = (v + \)uj. Thus, once 
(|B1|) is solved for x, we know the energy levels E v . 

There are two limiting cases corresponding to the har- 
monic potential (p = 1) and the half-harmonic potential 
(p = 0) that can be treated exactly. As 1/T(z) is an en- 
tire function in C which vanishes at all negative integers, 
we obtain: 

• For p = l,x = 3/4— n or x — 1/4— n, neff. Thus, 
i/ = n£N and E n = (n + ^)ui which is indeed the 
spectrum of the harmonic oscillator. 

• For p = 0, x = 3/4 - n, n € N. Thus, v = n + 1/4 
and E n — (n + which is also kown to be the 
spectrum of the harmonic oscillator on the half-line. 

In these two limiting cases, energy levels are exactly 
equidistant within the spectrum. The only noticeable 
change regards the Maslov index of the half- harmonic os- 
cillator which is 3/4 due to the presence of a "wall- type" 
singularity in x — (see for example .24]). 

However, for p > 0, no such singularity occurs. We 
expect the spectrum to be asymptotically given by the 
semiclassical energy levels, that is E n ~ (n + ^)uj as 
ri — * oo. As it has already been shown the spectrum 
of the split harmonic potential is not regularly spaced 
for < p < 1. Approximate expressions for the energy 
levels are given in Q when p ~ 1 or p ~ |43| . Explicit 
corrections are given for the first five energy levels when 
p ~ 1 and it is noticed that "The striking feature of these 
results is that they alternate in sign". 

To explain this interesting phenomenon, we perform an 
asymptotic expansion of the energy levels E n for p > 
as n — * oo. Starting from I|B1|) and after some algebra, 
we obtain 



E n 



1 



Xn 



where 

Xn 



2 (n+lY 

(l + p) 3 (l-p ) 
1287rp 2 

(l + p) 3 (l-p) 



lo , pn — ► oo 



(B2) 



cos 



2tt / 




n -\ 




U + P V 





128ttp 2 



COS 



27Tp 
1+P 



(B3) 



Clearly, (|B2p indicates that E n converges asymptotically 
to the harmonic levels for all p > 0. This could seem 
surprising since we know that they don't for p = (they 
are equal to (n+ |)w instead). In this respect, this limit 
is clearly singular. Hence, the requirement pn —> oo for 
(lB2l . 

The second remark is that, compared to Family I and 
Family II potentials whose levels converge like £T 5 / 2 to- 
wards the harmonic ones, the levels of the split-harmonic 
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FIG. 7: Parametric evolution of the ninth energy level, Eg, of 
a split-harmonic oscillator with frequency ratio p with respect 
to the variable f = (1 — p)/(l + p) (o> = 1). The solid line 
indicates the exact numerical result obtained from IBH and 
the dashed line indicates the asymptotic result IIB2II [to be 
compared with Figure 3 of reference 



FIG. 8: Corrections to the harmonic levels vs the energy level 
n > Xn = (n + ^) 2 [E n — (n+ §)], for a split-harmonic potential 
with frequency ratio p — 0.9. Cross symbols represent the 
exact numerical result obtained by solving <BH and circles, 
the asymptotic result \n given by iH3l . 



oscillator converge like E~ 2 instead. For isochronous po- 
tentials with large parameter [3 (that is, close to a split- 
harmonic oscillator), we expect a crossover between an 
asymptotic E~ 2 convergence law in the low/medium part 
of the spectrum and a E~ 5 / 2 convergence law in the large 
energy limit. 

Finally, eq. (|B3|) explains the oscillatory behaviour of 
the split-harmonic levels noticed by F.H Stillinger and 
D.K Stillinger in Q as both n and p vary since they are 
involved in the cosine function of l)B3jl . In particular, it 
can be checked that Figure 3 of [l| is reproduced perfectly 
well by (|B3|I (see Fig. And the n/2 oscillations of 
E n (p) within the range p € [0,1] are also explained by 
(1531 . 

In Figs. |Sland|!51 we give two typical corrections to the 
harmonic spectrum for p = 0.9 and p = 0.3 respectively. 
We define X^ xact = (rt +^ ) 2 [£'° xact — (n+|)] and compare 
it with Xn given by (|B3|) (u> = 1). Exact values for the 
energy levels are obtained by solving (|B1|I numerically. 
We see that x^ xact and x« are in good agreement even for 
small values of n. But as a general trend, this agreement 
degrades as p — * 0. 



Exact 
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FIG. 9: Same as Fig. 03 but for p = 0.3. 



where 



(C2) 



APPENDIX C: FROM h(E) TO S(X). 

In this last appendix, we investigate the inverse prob- 
lem which consists in recovering the potential (or at least 
the S function) corresponding to a prescribed function 
12(E). Starting from expression l|52|) for 12(E), we have 



Now, using the properties of Abel transforms (cf. for 
example |2(||3(j), we invert IjClfl to obtain 



12 d f E 2 I 2 (E) 

3(«) = -to-3- / 7i — ^m dE - ( C3 ) 



h 2 u>dv J (v-E) 1 / 2 ' 




(E-v) 1 / 2 



dv = — 



12^E 2 I 2 {E) 
h 2 oj 



(CI) Provided the integral on the r.h.s. of (|C3(I can be evalu- 
ated, we can use (|C21 to calculate S(X). 
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1. First example: 12(E) = est. 

Let us assume that h(E) = h = est. The r.h.s. 
3j yields 



of 



32/; 



a 2 < 



Reinstating this expression in 1C2J1 . we find 



S 2 



32/ 2 /v 2 
"fi 2 w V 2 



(C4) 



(C5) 



where a and b are two constants of integration. Since we 
are dealing with analytic potentials, 5 is continuous in 
and as it is odd, 5(0) = 0. Passing to the limit v — ► in 
(|C5J| shows that b = 0. Hence, 



v 5/2 -£j {diV^/uj)} ~ v~ 3 / 2 as v -> 00 and the Mellin 
transform M[#i,l] = f£%(v)dv ^ 0. Thus, according 
to formula (4.10.25) of 26], the first term of JUJ scales 
like ZT 7 / 2 as E -> 00. It can be verified in the same way 
that the second term scales like E~ 9 / 2 and is negligible 
when compared to the first one. 



1 - S 2 



S2v 



h 2 LO 1 

16/ 2 v + 2a 



(C6) 



Now, the condition 5(0) = determines the last con- 
stant, a = — h 2 u>/ (32I2) and 



5 



>2v 



i6/ 2 



(C7) 



which determines S(X) for X > and, as S(X) — 
—S(—X), S(X) on the entire real line. Finally, the re- 
quirement |5(X)| < 1 leads to the conclusion that I2 has 
to be negative. Let us write I2 = —h 2 /3 2 /(8uj) for some 
parameter j3 > 0, ilC7(l simplifies to 



S(X) = 



(3X 



^1 + (3 2 X 2 



(C8) 



which is nothing but the 5 function of an isotonic poten- 
tial with scaling parameter (3 (see l|2"2"jl. (l2^|) and JTSJl )■ 



2. Second example: J a (_E) = - ± {uj2 h + ^ E)9/2 ■ 

Repeating what has been done in the previous example 
for the function 12(E) chosen above yields the following 
5 function 

2X[35 + 42X 2 + 15X 4 ] 1/2 
1 ' [105 + 455X 2 + 483X 4 + 165X 5 } 1 / 2 ' 1 J 

This function is readily odd and VX € M, |5(X)| < 1. 
The corresponding potential has no singularity on the 
real line. 

For (|C9|) , the asymptotic behaviour of the fourth WKB 
term 14(E), given by J55j), is proportional to E~ 7 1 2 . To 
see this, we reinstate <|C9t into the function G\(X) de- 
fined in l|54|l and we perform an asymptotic analysis 
of the first term of l|55|) . It turns out that, g\(v) := 
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